Particle creation in an oscillating spherical cavity 
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Abstract 

We study the creation of massless scalar particles from the quantum vacuum due to the 
dynamical Casimir effect by spherical shell with oscillating radius. In the case of a small 
amplitude of the oscillation, to solve the infinite set of coupled differential equations for the 
instantaneous basis expansion coefficients we use the method based on the time-dependent 
perturbation theory of the quantum mechanics. To the first order of the amplitude we derive 
the expressions for the number of the created particles for both parametric resonance and 
non-resonance cases. 
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1 Introduction 



The Casimir effect is one of the most interesting manifestations of nontrivial properties of the 
vacuum state in a quantum field theory (for a review see [[j], ||, |3], [|) and can be viewed as 
a polarization of vacuum by boundary conditions. A new phenomenon, a quantum creation 
of particles (the dynamical Casimir effect) occurs when the geometry of the system varies in 
time. In two dimensional spacetime and for conformally invariant fields the problem with dy- 
namical boundaries can be mapped to the corresponding static problem and hence allows a 
complete study (see [||, Q and references therein). In higher dimensions the problem is much 
more complicated and only partial results are available. The vacuum stress induced by uniform 
acceleration of a perfectly reflecting plane is considered in |5|. The corresponding problem for 
a sphere expanding in the four-dimensional spacetime with constant acceleration is investigated 
by Frolov and Serebriany |7| in the perfectly reflecting case and by Frolov and Singh [|| for 
semi-transparent boundaries. For more general cases of motion the problem of particle and 
energy creation is considered on the base of various perturbation methods ||]-|23|](in the case of 
plane boundaries for more complete list of references see |2(J). It has been shown that a gradual 
accumulation of small changes in the quantum state of the field could result in a significant 
observable effect. A new application of the dynamical Casimir effect has recently appeared in 



connection with the suggestion by Schwinger [24| that the photon production associated with 



changes in the quantum electro dynamic vacuum state arising from a collapsing dielectric bubble 
could be relevant for sonoluminescence. For the further developments and discussions of this 
quantum-vacuum approach see p5f|-|3l[| and references therein. 



In our previous paper [32] we studied the problem of particle creation from the quantum 
vacuum by a spherical shell with time-dependent radius. We have considered examples for 
which the mean number of particles can be explicitly evaluated in the adiabatic approximation 
when the squeezing effect is dominant. In the present paper the case is considered when the 
sphere radius performs oscillations with a small amplitude and the expressions are derived for 
the number of created particles to the first order of the perturbation theory. 

2 Quantum scalar field inside a sphere 
with time-dependent radius 

Consider a scalar field <p satisfying Dirichlet boundary condition on the surface of a sphere with 
time-dependent radius a = a(t): 

( d 2 \ 

[ Qfi ~ A ) POM) = °> Plr=a(t)=°- (1) 

Following Ref.[32] we expand the corresponding eigenfunctions for the interior region in a series 
with respect to the instantaneous basis: 



Wmn{X) = \/ -37^ 2^ 5p \ Y lm (6,(p), (2) 

V a \ l ) j. = i Ji\Ji,k) 
I = 0,1,2,..., -l<m<l, n=l,2,..., (3) 

where (r, 8, cp) are standard spherical coordinates, ji k is the n-th zero for the spherical Bessel 
function ji(z), ji(ji n ) = 0) Yi m (9,ip) the spherical harmonic. These functions automatically 
satisfy boundary condition ([!]). Substituting (g) into field equation (p]) we arrive at an infinite 
set of coupled differential equations |32] 



Q l nk + ufk (tWnk = 2h J2 tfnp a pk + ^ iip a pk + h% J2 iip a P i a kj » ( 4 ) 
p=l p=l P,j=l 
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where overdot stands for the time derivative, 



h = -, u lk (t) = (5) 

and the time dependent antisymmetric coefficients a l nk are determined as 

l _ j 0, k = n 

ank ~\ -IjK.Mj/ijL .ih) k + n ■ W 

If the sphere is asymptotically static at past and future then the in- and out- vacuum states 
can be defined by using the solutions for coefficients corresponding to the in- and out-modes 
v!mn(*) with asymptotics 



ik\t) - ^r=S nk , t^-oo (7) 



e In 



gST(t) - ^=5n fc , f-+oo, (8) 

where we use the notations 

u>Z = ^, uC = —, o± = lim a(t) (9) 
for the corresponding eigenfrequencies. The expansion coefficients for the in- and out-modes are 



related by the transformation [32] 



where a l n j and • are Bogoliubov coefficients for the Bogoliubov transformation between the 

in- and out-modes V^mn^) and (£)■ In particular, taking into account (Q), in the limit 

t — > +00 from (|l0| ) we receive 



=={ a nk e lk + Pnk e lk ) 
, .out \ ' 



'ink — 

The mean number of out-particles produced in a given mode with quantum numbers l,m,n in 
the in- vacuum state is determined by the Bogoliubov coefficient [3 l nk : 

oo 

< in|iV /mn |in >= \Plk\ 2 , (12) 

k=l 

and the total number of created scalar particles is obtained by taking the sum over all the modes 

oo oo 

<in|iV|in>=£(2/ + l) £ \[3 l nk \ 2 . (13) 

1=0 n,k=l 

From the form of equations (||) it follows that there are two types of effects which lead to 
the particle creation (see also [19Q. The first one, called squeezing of the vacuum, is due to the 
nonstationary eigenfrequencies coi n (t) as a result of a dynamical change of the radius of the sphere 
and is described by the second term on the left of (|j). The second one, referred as acceleration 
effect, is due to the motion of the boundary and comes from the terms on the right of (|j). In [32] 
we have considered the particle creation inside a sphere in the adiabatic approximation when 
the squeezing effect is dominant. In the next section we will consider another approximation. 
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3 Particle creation by an oscillating sphere 

Below we will consider the case when the sphere radius performs small oscillations at a frequency 
f2 during a period T: 



a(t) 



ao[l + esin(nt)], 0<t<T 

ao, t < 0, t > T ' 



(14) 



.out 



(0) 
In 



where ao is the mean radius, e is a small parameter. For this type of motion ujf^ = u>Y"" = u 
jln/ a o- We will assume that for t < the field is in the in- vacuum state. This means that we 
need to solve Eq.(H) with the initial condition 



(o) , 

i & in 

Qnkit) = ^=S nk , t < 0. 



2u> 



(o) 

In 



(15) 



To derive the particle number created we will follow the scheme developed in [17] for the 
case of a perfect plane cavity with vibrating walls. First of all we note that introducing new 
functions 



x nk±(t) 



> u 



(0) 

Ik 



<L ■ ' '' 



(0) 
Ik ■ 



(16) 



to the first order of e the set of equations (|j) can be replaced by a coupled first-order differential 
equations 



X l nk± = ±iJ? k ) Xl k± ^iJ* ) e S \n(nt)(x l nk+ + Xl 



nk- 



(17) 



eft 



UJ 



(0) 
Ik 



cos(m) «5*v w 5 } " xl ^- 



± — , sin 



UJ 



(0) 
Ik 



3 k ( vl , vl 



UJ 



Introducing the column vector 



( X L- \ 



X n2- 



\ 



(18) 



this system can be written in the matrix form 



dt 



X l n (t) = V^X l n {t) + eV^X l n {t). 



Here the components for the matrices and are 

fA°) ■ (°) X X T/W V" 

where 

1 

J ik°kj 



isQt 



s=± 



V 



ka,ja' 



\ 



UJ 



(0) 



2 ' ^ 



(19) 



(20) 



(21) 
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and s, a, a' = +, — . To solve Eq. (|T^ ) we use a perturbation expansion 

where at the zeroth order, by taking into account initial condition (|l5l), one has 



X. 



(0)1 



• (o) , 

n hdrr-e 'k . 



v nka ~ °nkOa- 

The first order term can be easily found on the base of this expression: 



(22) 
(23) 



X 



nkcr 



■ (0). 



dt'e 



,(0), 



Vy(l) X W 
/ j ka,ja' nka' '" 



(0)1 



3-P 



(24) 



Using (|23|) and (|2C|) one finds 
where we use the notation 



(0) =F 



^Ik 
.(0) 



(0) 



(25) 



(0) 

E2 k (t) 



(i/m) 



■ (o). 
-i(m+o),v )t 



■ (o) 4 



for m = 
for m/0 



(26) 



For the instantaneous basis expansion coefficients this yields 



■ (o) 4 

-tU)} u 't 



Q l nk(t) 



2lo 



2lo 



(0) 
ik 



k— ,ra- 



(0) 



(0) 



,/ 0) 



J k-,n-^ .,(0) o ,.,(0) 



y fe-,n- 



J k-,n- 



,, (0) 
TP ^ 
A ,(0). 



i fe in 



+ 0(e 2 



(27) 



This expression includes terms proportional to t which are due to the parametric resonance. In 
the situation ujf^t » 1 the resonance terms are dominant and solution (|27|) becomes 



Qnk(t) 



■ (o). 



2w 



+ v 



fc — ,71- 



■ (0). 



2w 



(o) 
Ik 



+ jo;,, 'tr 



,(0) 



(28) 



(0) o , (0) + V k- 
Zfc ' In 



■ (.0). 



ik l s 



(0) 



(()) 



in ' Ik 



After the time interval T the sphere is static with radius ao and the function q l n f.{t) has the form 
(11). From the continuity condition between the functions (11) and (28) at t = T one has 



a 



nk 



0; 



'nk 



Kk + £T 

-eTv' 



V h _ n _fi (0) . (0) + Vr_ 5 (0) „ . (0) 
Ik ' In In ' Ik 



t k ' in 



By taking into account expressions (^l|) and @ one has 

1 



cT ./, ,(0) (0), 



(0). 

/ 1 J 



(29) 
(30) 

(31) 



Now from Eq.(12) for the number of out-particles with the energy and quantum numbers 



I, m created by the parametric resonance in the in- vacuum state we receive 



1 



< in|iV Zmn |in >= -{eTf ^ u\ n >u\ k } 5 {0) (0) . 

t-t , _, Ik ' In 

k=l 



(32) 
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As we see the necessary condition for the parametric resonance is 



n 



— (jl,q + jl,p) 

a 



(33) 



for some values q and p. 

Now we turn to calculating the Bogoliubov coefficient (3 l nk in the non-resonance case. By 
using the expressions (|26| ) from the continuity condition between the functions ( |Tl"l) and 
and their derivatives at t = T after some algebra we get 



A 



'nk 



2,exp - l T{^+^),2\ £ 



(0) 



JaflT/2_ 



sin 



For the number of quanta we derive from here 



(0) 

In 



+ UJ, 



(0) 

Ik 



an 



(34) 



< in|iV; mri |in > 



k=l 



'k—,n— 



2 sin" 



(0) 



{ U ln 



(0) , ,,(o) 



+ W, 



/A- 



aft) 



+ 



V k-,n- V k-,n- C0S(CIT) 



(J, 



(0) 



n 2 



COS(OT) — COS T ( Ufa + LU 



(35) 



Note that in the resonance case from the term with a = + of Eq.(34) we recover the result (|30[) . 

Up to now we have considered the scalar field with the Dirichlet boundary condition. In 
a similar way it is easy to generalize the corresponding results for other type of boundary 
conditions. For example, in the case of the Neumann scalar the number of particles created 

n = 1, 2, ... are 



inside an spherical shell is given by expressions ([32;) and (35) where now au\ n 



,(0) 



zeros for the function j[{z), j[{auj\^) = 0. 



4 Concluding remarks 

Our interest in this paper has been in the quantum radiation by an spherical mirror with time- 
dependent radius given by (jlj) and under assumption that the amplitude of oscillation is small. 
To define and interpret asymptotic spaces of physical states uniquely we assumed that the 
oscillation of the sphere lasts a finite period of time. We consider a massless scalar field obeying 
the Dirichlet boundary condition on the sphere surface. To solve the corresponding equations for 
the spontaneous basis expansion coefficients in the short-time limit, eQT <C 1, we use the method 



developed in [17] for the case of plane boundaries. To the first order in perturbation theory we 
derive the expression for the number of particles created inside the sphere in both resonance and 
non-resonance cases. In the resonance case the oscillation frequency is given by expression (|33|). 
The corresponding solution includes resonance terms proportional to t and oscillating parts. In 
the situation, when QT 3> 1 the resonance terms are dominant, and the total number of particles 
created in the mode with quantum numbers Iran and with frequency oJ^ is given by (|32|). In 
this case the number of particles grows quadratically in time. In the non-resonance case the 
corresponding formula for the number of particles has the form (|35|). The generalization of these 
formulas for the other boundary conditions is straightforward. Another possible generalization 
is to perform the higher order calculations with respect to e. We hope that the results obtained 
in present paper may have an application in studies of quantum radiation created by bubble 
formation during first order phase transitions in the Early Universe. 



6 



References 

[1] G. Plunien, B. Mueller, and W. Greiner, Phys. Rep. 134, 87 (1986). 

[2] V. M. Mostepanenko and N. N. Trunov, The Casimir effect and its applications (Oxford 
Science Publications, New York, 1997). 

[3] K. A. Milton, in Applied Field Theory, ed. C. Lee, H. Min, and Q-H. Park (Chungbum, 



Seul, 1999) p.l, |hep-th/9901011 



[4] N. D. Birrell and P. C. W. Davies, Quantum fields in curved space (Chambridge University 
Press, 1982). 

[5] P. Candelas and D. Deutsch, Proc. Roy. Soc. A354, 79 (1977). 

[6] V. P. Frolov and E. M. Serebriany, J. Phys. A: Math. Gen. 12, 2415 (1979). 

[7] V. P. Frolov and E. M. Serebriany, J. Phys. A: Math. Gen. 13, 3205 (1980). 

[8] V. Frolov and D. Singh, Class. Quantum Grav., 16, 3693 (1999). 

[9] L. H. Ford and A. Vilenkin, Phys. Rev. D 25, 2569 (1982). 

[10] G. Calucci, J. Phys. A: Math. Gen. 25, 3873 (1992). 

[11] G. Barton and C. Eberlein, Ann. Phys.,N.Y., 227, 222 (1993). 

[12] C. K. Law, Phys. Rev. A 49, 433 (1994). 

[13] R. Jauregui, C. Villarreal, and S. Hacyan, Mod. Phys. Lett A10,7 (1995). 

[14] E. Sassaroli, Y. N. Srivastava, and A. Widom, Phys. Rev. A 50, 1027 (1994). 

[15] V. V. Dodonov and A. B. Klimov, Phys. Rev. A 53, 2664 (1996). 

[16] A. Lambrecht, M.-T. Jaekel, and S. Reynaud, Phys. Rev. Lett. 77, 615 (1996). 

[17] Jeong- Young Ji, Hyun-Hee Jung, Jong-Woong Park, and Kwang-Sup Soh, Phys. Rev. A56, 
4440 (1997). 

[18] L. Sh. Grigoryan and A. A. Saharian, Izv. Akad. Nauk Arm. SSR Fiz. 32, 223 (1997); 32, 
275 (1997) (Sov. J. Contemp. Phys.). 

[19] R. Schutzhold, G. Plunien, and G. Soff, Phys. Rev. A 57, 2311 (1998). 

[20] V. V. Dodonov, J. Phys. A: Math. Gen. 31, 9835 (1998). 

[21] R. Golestanian and M. Kardar, Phys. Rev. A 58, 1713 (1998). 

[22] L. Hadasz, M. Sadzikowski, and P. Wegrzin, "Quantum Radiation from Spherical Mirrors", 
|hep-th/9803032| . 

[23] W. G. Andersen and W. Israel, Phys. Rev. D 60, 084003 (1999). 

[24] J. Schwinger, Proc. Nat. Acad. Sci. 90, 985, 2105, 4505, 7285 (1993); 91, 6473 (1994). 

[25] C. Eberlein, Phys. Rev. A53, 2772 (1996); Phys. Rev. Lett. 76, 3842 (1996). 



7 



[26] K. A. Milton, "Casimir Energy for a Spherical Cavity in a Dielectric:Toward a Model for 
Sonoluminescence" , tiep-th/9510091. 



[27] C. E. Carlson, C. Molina-Paris, J. Perez-Mercader, and M. Visser, Phys. Rev. D 56, 1262 
(1997). 

[28] C. Molina-Pari's and M. Visser, Phys. Rev. D 56, 6629 (1997). 
[29] K. Milton and J. Ng, Phys. Rev. E 57, 5504 (1998). 

[30] S. Liberati, M. Visser, F. Belgiorno, and D. W. Sciama, J. Phys. A: Math. Gen. 33, 2251 
(2000). 

[31] S. Liberati, M. Visser, F. Belgiorno, and D. W. Sciama, Phys. Rev. D 61, 085023 (2000). 

[32] M. R. Setare and A. A. Saharian, "Particle creation by moving spherical shell in the dy- 
namical Casimir effect", hep-th/0101149 ; to appear in Mod. Phys. Lett. A. 



8 



